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Abstract 

Wind tunnel measurements of the rotor trim, blade airloads, and structural loads of a full-scale UH-60A 
Black Hawk main rotor are compared with calculations obtained using the comprehensive rotorcraft analysis 
CAMRAD II and a coupled CAMRAD II/OVERFLOW 2 analysis. A speed sweep at constant lift up to 
an advance ratio of 0.4 and a thrust sweep at constant speed into deep stall are investigated. The coupled 
analysis shows significant improvement over comprehensive analysis. Normal force phase is better captured 
and pitching moment magnitudes are better predicted including the magnitude and phase of the two stall 
events in the fourth quadrant at the deeply stalled condition. Structural loads are, in general, improved with 
the coupled analysis, but the magnitude of chord bending moment is still significantly underpredicted. As 
there are three modes around 4 and 5/rev frequencies, the structural responses to the 5/rev airloads due to 
dynamic stall are magnified and thus care must be taken in the analysis of the deeply stalled condition. 


Nomenclature 

Cl rotor lift coefficient 

Ct rotor thrust coefficient 

M 2 Cn nondimensional section normal force 

M 2 Cm nondimensional section pitching moment 

a c wall-effect corrected shaft angle 

a s geometric shaft angle (positive for rearward tilt) 

/ 3„ coning, deg 

Pic longitudinal flapping, deg 

Pis lateral flapping, deg 

H advance ratio 

a solidity, 0.0826 

d 0 collective, deg 

die lateral cyclic, deg 

d\ s longitudinal cyclic, deg 


Introduction 

A full-scale wind tunnel test of a UH-60A rotor was recently 
conducted (May 2010) by NASA and the U.S. Army in the 
National Full-Scale Aerodynamics Complex (NFAC) 40- by 
80-Foot Wind Tunnel to provide an expanded database for 
validating and enhancing rotorcraft aeromechanics prediction 
methodology [1], This test acquired a comprehensive set 
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of validation-quality measurements on a full-scale pressure- 
instrumented rotor system at conditions that challenge the 
most sophisticated modeling and simulation tools. Figure 1 
shows these UH-60A rotor blades installed on the NFAC 
Farge Rotor Test Apparatus (FRTA) in the wind tunnel test 
section. 

Data were acquired over a wide range of operating conditions, 
including conditions of airspeed up to 175 kt, thrust up to 
32,000 lb, and advance ratio up to 1 .0. The database provides 
aerodynamic pressures, structural loads, control positions, 
and rotor forces and moments, allowing for the validation of 
both aerodynamic and structural models. These unique and 
extensive test data sets provide a useful resource that can be 
used to examine the rotor behavior in the most challenging 
conditions. 

For the present investigation, two parametric sweeps were 
selected for analysis. First, a speed sweep was selected 
to test simulation accuracy over a wide range of advance 
ratios. Second, a thrust sweep was selected to test simulation 
accuracy under conditions extending from unstalled through 
some deep stall conditions. The rotor blade aerodynamic 
environment at high speed is characterized by compressibility 
and a poor distribution of loading over the rotor disk. At 
high thrust, aerodynamics are dominated by dynamic stall and 
large negative pitching moment. 

Romander et al. [2] provided an initial correlation of airloads 
with measured data for the same test sets considered in 
the present paper using the coupled CAMRAD II [3] and 
OVERFFOW 2 [4], In the present paper, both CAMRAD II 


1 



alone and coupled CAMRAD II and OVERFLOW 2 analyses 
are carried out and rotor structural loads as well as airloads 
results are compared with the UH-60A rotor wind tunnel test 
data. Detailed examination and parametric studies are also 
conducted for the deeply stalled condition, where the previous 
analysis [2] had difficulty getting good correlation. 

Description of the Test 

The experiment was conducted in the NFAC 40- by 80-Foot 
Wind Tunnel using a UH-60A rotor system mounted on the 
NFAC LRTA. The blades used in this test program were 
the same 4 rotor blades flown during the UH-60A Airloads 
Program [5], A detailed description of the experiment can be 
found in Ref. 1, including information on the test hardware, 
instrumentation, data acquisition and reduction systems, rotor 
control systems, and standard test procedures. Information 
relevant to the current study is provided below. 

Instrumentation 

Two of these blades were heavily instrumented: one with 242 
pressure transducers and one with a mix of strain-gages and 
accelerometers. The pressure blade was originally built with 
242 sub-miniature pipette type absolute pressure transducers 
embedded below the skin surface of the blade. Absolute 
pressures were measured at nine radial locations (r/R = 0.225, 
0.40, 0.55, 0.675, 0.775, 0.865, 0.92, 0.965, and 0.99) as 
shown in Fig. 2(a). Blade section normal force, pitching 
moment, and chord force were obtained by integrating the 
absolute pressures. For the blade airloads presented in 
this paper, only working transducers were included in the 
integrations. Detailed examination of pressures at each 
transducer are still in progress to ensure the accuracy of the 
integrated airloads. 

Blade flap bending, chord bending, and torsion moments were 
measured with two- or four-leg strain-gage bridges bonded 
to the second instrumented blade. The gages were located 
at the blade root (11.3% radius) and then evenly distributed 
along the blade at 10% increments of the rotor radius (20% to 
90%) as shown in Fig. 2(b). Most of the gauges were bonded 
directly to the blade spar and used during the Airloads flight 
test. Five additional torsion gages were bonded to the blade 
skin, including one at the same radial station as a spar-bonded 
gage to ensure measurement consistency. 

The blade root motions occur around elastomeric bearings 
and the “hinges” are the focal points of the bearings. Two 
independent, specially designed blade motion measurement 
systems were used to determine blade flap, lag, and pitch 
angles at the root of each blade. The first system, 
designated the Blade Motion Hardware (BMH) or “crab 
arm”, was used during the Airloads flight test. The BMH 
are mechanical measurement devices that span the main 
rotor blade hinge points and include three Rotary Variable 


Differential Transformers (RVDTs) on each blade to measure 
the relative motions of its own articulations. The true blade 
motions are obtained through three kinematic equations that 
account for the cross-coupling between measurements. The 
second system, designated the Rotor State Measurement 
System, is composed of four sets of three laser distance 
transducers (one set mounted to each hub arm). Each 
transducer produces a laser beam that is reflected back to 
a sensor in the same enclosure by reflective objects within 
a specified distance range. The transducer then produces a 
voltage proportional to the distance to the object (and related 
to one of the blade root angles). The accuracy of these blade 
root motion measurement systems is currently under review. 

Test procedure 

A level flight speed sweep was conducted to provide data 
that showed the effects of advance ratio for representative 
flight conditions. For each wind tunnel test condition, rotor 
collective, cyclic, and shaft angle were adjusted to match 
the target lift, propulsive force, and hub moments. The 1-g 
level flight speed sweeps were performed at three lift levels, 
Cl/ct = 0.08, 0.09, and 0.10, up to an advance ratio of 0.4. 

A thrust sweep was performed using specified rotor thrust and 
zero hub pitch and roll moments as trim targets until stall 
was detected. Then, collective pitch was manually varied to 
better define the phenomena. During the collective sweeps 
only zero hub moments were achieved (thrust not trimmed). 
The thrust sweeps were conducted at three different tip Mach 
numbers; Mu p = 0.625, 0.650, and 0.675. The majority of 
the data were acquired at the baseline M tlp = 0.650 with 
various shaft angles and advance ratios. A limited number 
of sweeps at M t i P = 0.625 were conducted to attain higher 
non-dimensional thrusts and advance ratios without reaching 
load limits. 

In this paper, the speed sweep ranged from /i = 0.15 to 0.4 
with a constant Cl /er = 0.09 and a tip Mach number of 0.65 
(Run 52). The selected thrust sweep was at // = 0.3 with a tip 
Mach number of 0.625 (Run 45). The test points are plotted 
in Fig. 3 along with aerodynamic rotor lift boundary obtained 
in wind tunnel testing of a model rotor by McHugh [6], The 
maximum thrust obtained in this test is slightly higher than the 
McHugh lift boundary. Among those sweeps, five data points 
(circled in Fig. 3) were further selected for detailed analysis. 
From the speed sweep, three points were selected; // = 0.2 
(Run 52, Point 20), // = 0.3 (Run 52, Point 31), and // = 0.4 
(Run 52, Point 52) at Cl/ct = 0.09. From the collective 
sweep, two points including the extreme thrust condition were 
selected; Ct / cr = 0.06 (Run 45, Point 30) and Ct /cr = 0. 1255 
(Run 45, Point 38) at /i = 0.3. 
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Description of Analytical Method 

The analytical results were obtained using the comprehensive 
analysis CAMRAD II alone and separately using coupled 
CAMRAD II/OVERFLOW 2. This section describes each 
method and how they are coupled to produce a higher fidelity 
solution. 

CAMRAD II 

The analytical results were obtained using the comprehensive 
analysis CAMRAD II. CAMRAD II is an aeromechanics 
analysis of rotorcraft that incorporates a combination 
of advanced technologies including multibody dynamics, 
nonlinear finite elements, and rotorcraft aerodynamics [3]. 
CAMRAD II has been used extensively for correlation of 
performance and loads measurements of the UH-60A in 
various flight conditions [7-10], The aerodynamic model 
is based on second-order lifting line theory [11]. The 
blade section aerodynamic modeling in lifting line theory 
is unsteady, compressible, viscous flow about an infinite 
wing in a uniform flow consisting of a yawed freestream 
and wake-induced velocity. This problem is modeled within 
CAMRAD II as two-dimensional, steady, compressible, 
viscous flow (airfoil tables), plus corrections for swept and 
yawed flow, spanwise drag, unsteady loads, and dynamic stall. 
The wake modeling of lifting line theory is an incompressible 
vortex wake behind the lifting line with distorted geometry 
and rollup. The wake analysis calculates the rotor nonuniform 
induced velocities using rigid, prescribed, or free wake 
geometry. Two wake analyses are available in CAMRAD II: 
a rolled-up model and a multiple-trailer (with consolidation) 
model. The rolled-up wake model is based on the assumption 
that a tip vortex forms at the outboard blade tip. Because 
of its simplicity and efficiency, the rolled-up model has long 
been used for helicopter rotors . The multiple-trailer mode 1 has 
also been available and, with the consolidation feature, has 
been applied recently with success to tiltrotor and helicopter 
airloads calculations [9, 12]. 

The multiple-trailer model has the far wake trailed vorticity 
divided into several spanwise panels to provide more detailed 
structure for the inboard vorticity. This model has a discrete 
trailed vortex line emanating from each of the aerodynamic 
panel edges. The calculation of the free wake geometry 
includes the distortion of all of these trailed lines. With 
multiple far wake trailed vorticity panels, the trailed lines at 
the aerodynamic panel edges can be consolidated into rolled- 
up lines using the trailed vorticity moment to scale the rate 
of roll-up. It is assumed that all the vorticity in each set 
eventually rolls up into a single vortex located at the centroid 
of the original vorticity distribution. 

In this work, an isolated rotor was modeled as a flexible 
blade with 11 nonlinear finite elements and 21 aerodynamic 
panels. Detailed rotor control system geometry, stiffness. 


and lag damper were also incorporated. When not coupling 
to OVERFLOW, a dual-peak rigid wake model was used, 
except at /i = 0.2 where the free-wake model was used. The 
multiple-trailer with consolidation model was also used for 
the = 0.2 case and the results were compared with those 
with the baseline rolled-up wake model. 

OVERFLOW 2 

All Navier-Stokes Computational Fluid Dynamics 
(CFD) analyses presented herein were performed using 
OVERFLOW 2 version 2.2b [4], OVERFLOW 2 is an 
overset, structured-mesh flow solver developed at NASA. 
For two decades the OVERFLOW solver has served to 
analyze a variety of rotorcraft under a wide range of flight 
conditions [13]. OVERFLOW 2 offers a wide variety 
of numerical schemes, turbulence models, and boundary 
conditions. For the present study, OVERFLOW 2 was run 
with 4th order central differencing and 4th order artificial 
dissipation in space. Time marching was performed using 
a 2nd order dual timestepping scheme. Turbulence was 
modeled near blade surfaces using the Spalart-Almaras one 
equation model with rotational corrections. The turbulence 
model was deactivated in regions one chord length or 
further from the rotor blades to reduce numerical dissipation 
of the wake. Blade surfaces were modeled as viscous, 
adiabatic walls; outer boundaries were modeled using a 
characteristic condition imposing freestream quantities. 
OVERFLOW 2 computes the flowfield by discretizing the 
Navier-Stokes equations on a series of overset, structured 
grids. Grids modeling the rotor blades were body-fitted 
and curvilinear. These grids, often called near-body grids, 
extended approximately one chord length from the blade 
surface. The near-body grids were nested within one or more 
grids, called off-body grids, which filled the space between 
the rotor and the boundary of the computational domain. The 
OVERFLOW 2 model included a notional hub, but the LRTA 
and wind tunnel struts were not modeled. All grids exchanged 
flow information in regions of overlap. The amount of this 
overlap was sufficient to support full 4th order accuracy at the 
boundaries. 

The near-body grid representing the bulk of each rotor blade 
had an ”0” topology with 157 points wrapping around the 
blade chordwise, 163 points along the blade span, and 75 
points normal to the surface. The initial spacing at the 
blade surface had a y + value of 1. This grid system is 
illustrated in Fig. 4(a). The near-body grid system (the hub 
and all four blades) contained approximately 11.6 million 
points. The off-body grid system used a series of ever larger 
Cartesian grids to create shells expanding outward from the 
near-body grid set. The grid point spacing within each shell 
is twice that of the shell immediately preceding it. The 
finest off-body grid had a spacing equivalent to 0. 1 tip-chord 
lengths in all three directions. Seven such shells created a 
cubic computational domain spanning ten rotor radii in every 
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direction. Domain boundaries in this scenario were set to a 
freestream characteristic condition thereby simulating a rotor 
operating in free air. Flow in the off-body grids was treated 
as inviscid and the turbulence model was deactivated. This 
off-body grid set consisted of 15.5 million grid points and is 
depicted in Fig. 4(b). 

The near-body and off-body grid combination totals 27.1 
million grid points. To reduce computation time, CFD 
simulations were run using a hybrid distributed/shared 
memory scheme with 160 Message Passing Interface 
(distributed memory) ranks and eight OpenMP threads 
(shared memory) per MP1 rank for a total of 1280 parallel 
tasks. OVERFLOW 2 required approximately 31 minutes 
to advance the solution for this configuration 1/4 of a rotor 
revolution using 1280 CPUs of an SGI Altix ICE computer. 

Coupled Analysis of CAMRAD II/OVERFLOW 2 

CAMRAD II uses a lower-fidelity aerodynamics model than 
that available in modern CFD codes, and most CFD codes 
lack the sophisticated Computational Structural Dynamics 
(CSD) and trim capabilities of comprehensive codes like 
CAMRAD II. Coupling a CFD code (e.g. OVERFLOW 2) 
to a comprehensive code (e.g. CAMRAD II) marries the 
strengths of the two approaches and produces the highest- 
fidelity solution currently possible. For this study, coupling 
is achieved by alternate execution of OVERFLOW 2 and 
CAMRAD II. At the end of each code’s turn to execute, 
it passes data to the next code. The data passed from 
OVERFLOW 2 to CAMRAD II is airload data integrated from 
its Navier-Stokes model of the UH-60 rotor. This airload 
data is used to augment CAMRAD II’s internal aerodynamics 
model (which consists of airfoil tables and a lower-order 
wake model). At the end of its execution, CAMRAD II 
generates updated control positions and a description of 
how the blade deforms elastically as it revolves around the 
shaft. These quantities are used to give OVERFLOW 2’s 
grids a realistic motion in response to the aerodynamic 
environment. This algorithm, called the delta coupling 
technique, was pioneered by Tung et al. [14] and implemented 
in OVERFLOW by Potsdam et al. [15] Significantly improved 
airloads prediction capability has been demonstrated for the 
UH-60A rotor in steady level flight conditions using this 
loosely coupled approach [2, 15, 16]. The CFD solution is 
advanced 1/4 revolution during a coupling iteration because 
this allows each of the rotor’s four blades to sweep through 
a full quadrant of the rotor disk. Taken in aggregate, the 
four blades thereby determine the airloads at every azimuth 
for every coupling iteration. Convergence of the coupling 
process was determined by monitoring blade airloads for 
periodicity. When the airloads did not vary significantly from 
one coupling iteration to the next, the solution was judged 
to be converged. For the present analysis, this generally 
occurred after 24 coupling iterations. Since OVERFLOW 2 
was allowed to iterate for 1/4 revolution between coupling 


exchanges, this equates to 6 full revolutions for the converged 
solution. A fully converged coupled solution required 
approximately 17 hours to compute on 1280 SGI Altix ICE 
processors. 

For coupled analysis, CAMRAD II input was identical 
to uncoupled CAMRAD II input except the wake model 
was set to uniform inflow. The airload data passed to 
CAMRAD II accounted for both sectional aerodynamics 
and wake effects. The lifting-line aerodynamics and wake 
model were essentially replaced by airloads provided by 
OVERFLOW 2. 

Results and Discussion 

In this section, selected data from the wind tunnel test are 
compared with predictions from the two analysis methods. 
These test data include rotor trim, blade hinge motion, normal 
force, pitching moment, flap bending, chord bending, and 
torsion moments for both speed and thrust sweeps. 

Trim and blade motion for speed sweep 

Figures 5 and 6 show the calculated and measured blade pitch 
and flapping angles at the hinge at three different advance 
ratios; // = 0.2 (Run 52, Point 20), /r = 0.3 (Run 52, Point 31), 
and fi = 0.4 (Run 52, Point 52) for the same C [J rr = 0.09. The 
trim solution for the speed sweep conditions in the analyses 
solves for the collective and cyclic angles required to match 
the measured thrust and hub moments with specified corrected 
shaft angle ( a c ). The measured hinge motions were obtained 
from the BMH hardware and the values averaged over each 
of the four blades were used for comparison. The measured 
collective angles increase as the advance ratio increases and 
the analyses show the same trends but overpredict by about 
1 deg at all speeds investigated. Note that different scales are 
used for each plot. The calculations show good correlation of 
cyclic pitch angles, except lateral cyclic angle at fi = 0.4 where 
the analysis underpredicted by about 2 deg. The coupled 
analysis shows slightly better agreement with the measured 
cyclic angles compared to CAMRAD II alone. 

The measured coning angles at the three different advance 
ratios are almost identical because lift is same for those 
conditions. The calculations show good agreement with 
the test data. The difference between the calculations and 
measured values are less than 0.2 deg. The difference between 
the two analyses is also small. 

Airloads for speed sweep 

Blade section normal force and pitching moment for the 
speed sweep are investigated in this section. Figure 7 shows 
nondimensional normal force at 92% radial station at three 
different advance ratios. At /r = 0.2, the measured normal 
force shows pulses due to blade vortex interaction (BVI) on 
both the advancing side and retreating side. The CAMRAD II 
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analysis significantly overpredicts the magnitude of the down- 
up impulse on the advancing side. The same overprediction 
was also observed with the UH-60A flight test correlation 
at [i = 0.15 [9]. This will be revisited in a later section. 
The coupled analysis accurately captures the magnitude of 
these pulses. At higher speed, the airloads on the blade 
tip region are generally characterized by negative lift at the 
end of the first quadrant and the beginning of the second 
quadrant [18] and both /j = 0.3 and 0.4 conditions show 
negative normal force at this radial location. This negative 
loading gets more severe as the advance ratio increases. The 
magnitude of the normal force is reasonably well predicted 
by CAMRAD II, but the phase angle correlation is fair to 
poor. The phase of the airloads from the coupled analysis is 
significantly improved over the CAMRAD II analysis. Small, 
high frequency oscillations in the test data in the first quadrant 
resulting from the wake interaction are also beginning to 
be captured in the coupled solution. Figure 8 shows 
nondimensional pitching moment at 92% radial station for the 
speed sweep. The CAMRAD II analysis shows a significant 
underprediction of advancing blade pitching moments at all 
advance ratios. Although there is an underprediction of 
peak-to-peak magnitude and an exaggerated dynamic stall 
cycle in the fourth quadrant at //, = 0.4, the coupled analysis 
significantly improves the overall correlation compared to the 
comprehensive code alone, especially the moment pulse due 
to BVI is better captured. 

Structural loads for speed sweep 

Figures 9 through 12 compare the calculated and measured 
structural loads for the speed sweep. Figure 9 shows 
oscillatory flap bending moment at 40% radial station. Steady 
values were removed from both test data and analyses. The 
CAMRAD II analysis showed poor correlation at n = 0.2 due 
to the overprediction of the magnitude of the down-up normal 
force impulse on the advancing side as shown in Fig. 7(a). 
The CAMRAD II calculated flap bending moments show fair 
to good correlation on magnitude at // = 0.3 and 0.4, but 
the phase differs significantly from the measurements. The 
coupled analysis shows better prediction capability than the 
CAMRAD II analysis. Waveform and phase correlation is 
improved in the first and second quadrants by the coupled 
analysis. 

Figure 10 shows oscillatory chord bending moment at 40% 
radial station. The measurements show significant high 
frequency content unlike flap bending moments. Both 
analyses underpredict the magnitude. Although the coupled 
analysis shows much sUonger high frequency content than the 
comprehensive analysis, especially at /r = 0.2 and 0.3, the 
peak-to-peak magnitude is significantly underpredicted and 
phase is not well captured. 

Figure 11 shows oscillatory torsion moment at 40% radial 
station. Correlation is poor at /j = 0.2 for both analyses. 


However both analyses, in general, show reasonably good 
correlation at higher advance ratios although peak-to-peak 
magnitude is underpredicted. The coupled analysis improves 
phase in the first and second quadrants and waveform in the 
third and fourth quadrants. The sharp drop around 250-deg 
azimuth is not well predicted even with the coupled analysis. 

Figure 12 shows oscillatory pitch link loads. The measured 
pitch link loads show similar waveform as the torsion 
moments at r/R = 0.4 on the advancing side, but much 
larger variation on the reUeating side. The magnitude is well 
predicted in the first and second quadrants by both analyses, 
especially at /i = 0.3 and 0.4, but the correlation is poor 
in the third and fourth quadrants. Although detailed shapes 
are better predicted by the coupled analysis, the sharp drop 
around 250-deg azimuth is not captured even with the coupled 
analysis. 

Figure 13 compares harmonic magnitude of flap bending, 
chord bending, and torsion moments along the blade span 
for /i = 0.3. Note that different scales are used for each 
harmonic. The coupled analysis shows better agreement with 
the measured data on the 2/rev flap bending moment and 
most of the 4 and 5/rev moments. Neither analysis does 
well for the 1, 2, and 3/rev chord bending moment. Even 
though a nonlinear lag damper model was used, both analyses 
significantly underpredict the magnitude of chord bending 
moment. 

Trim and blade motion for thrust sweep 

Figures 14 and 15 show the calculated and measured blade 
pitch and flapping angles at the hinge at two different thrust 
levels; C t /ct = 0.06 (Run 45, Point 30) and C T /cr = 0.1255 
(Run 45, Point 38) for the same //, = 0.3. The trim solution for 
the CAMRAD II was solved for the controls that produced 
rotor thrust and hub moments to match the measured values, 
with the rotor shaft angle of attack fixed at the measured 
values with correction. It should be noted that the measured 
hub moments are very small as zero hub moments were 
sought during the test. For the coupled analysis, a different 
trim strategy was used. First, baseline case was selected at 
Ct/ct = 0.08 and the simulation was trimmed to match the 
measured thrust. Collective deltas were derived from test 
data relative to this baseline and then applied to the simulated 
baseline to produce the remaining target points. This resulted 
in Ct/c = 0.0591 and 0.1243, respectively. For the maximum 
thrust condition, the estimated thrust is about 0.96% lower 
than the measured value, with collective of 13.34 deg. 

For the Ct / cr = 0.06 case, the coupled analysis shows 
slightly better collective and coning angle correlation with 
the measured data than the comprehensive analysis, but 
both analyses overpredict them. For the highest thrust 
condition, although collective is still about 1.5 deg higher 
than the measured value, the coupled analysis shows excellent 
agreement of coning angle. 
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Airloads for thrust sweep 

Figure 16 shows nondimensional section normal force at 
92% radial station for thrust sweep. The measured data 
show negative normal force at the end of first quadrant and 
beginning of second quadrant for both thrust levels. The 
measured normal force at Cy/cr = 0.06 shows BVI-induced 
pulses in the first quadrant, which is similar to that observed 
at lower advance ratio (Fig. 7(a)). At higher thrust, the 
measured data show two stall events in the fourth quadrant. 
At Ct/c = 0.06 condition, the coupled analysis captures 
the BVI-induced normal force variation in the first quadrant 
and significantly improves phase correlation on the advancing 
side. At Ct/cf = 0.1255 condition, the coupled analysis 
shows good correlation on the retreating side, capturing the 
large normal force fluctuations in the fourth quadrant due to 
dynamic stall events. However, the coupled analysis shows 
a large normal force variation in the first quadrant and this 
(along with pitching moment) has an important influence on 
high frequency structural loads that will be discussed later. 

Figure 17 shows nondimensional section pitching moment at 
92% radial station for the thrust sweep. The CAMRAD II 
analysis again shows a significant underprediction of 
advancing blade pitching moments. It should be noted that the 
CAMRAD II analysis did not use a dynamic stall model and 
thus there is no high frequency pitching moment variation at 
Ct/ct = 0.1255. The coupled analysis significantly improved 
the correlation by capturing BVI-induced high frequency 
impulse in the first quadrant at CtI<7 = 0.06 and the magnitude 
and phase of the two stall events in the fourth quadrant at 
the deeply stalled condition (Ct/ct = 0.1255). Similar to the 
normal force, the coupled analysis shows a pitching moment 
pulse in the first quadrant, which was not observed in the 
measured data. 

Structural loads for thrust sweep 

Figures 18 through 21 show structural loads correlation for 
the thrust sweep. The coupled analysis shows much better 
flap bending moment correlation, as shown in Fig. 18, by 
improving both magnitude and phase correlation with the 
measured data. For the chord bending moment shown 
in Fig. 19, both analyses underpredict the magnitude at 
Ct/c = 0.06, although the coupled analysis shows slightly 
better peak-to-peak magnitude. However, the coupled 
analysis shows significant 5/rev chord bending moment. This 
overprediction of the 5/rev harmonic is also observed in 
the torsion moment and pitch link loads shown in Figs. 20 
and 21, respectively. This can be explained from a coupled 
aerodynamic/structural dynamics point of view. Normal force 
and pitching moment pulses in the first quadrant, combined 
with two dynamic stall cycles in the fourth quadrant, generate 
a strong 5/rev forcing function. As the blade 2nd chord 
frequency is about 4.69/rev and first torsion frequency is about 
4.53/rev, the structural responses to the 5/rev airloads are 
magnified. This will be further examined in a later section. 


Effect of wake modeling at low speed 

The CAMRAD II analysis significantly overpredicted the 
magnitude of the down-up impulse on the advancing side at 
/r = 0.2 as shown in Fig. 7(a). A similar trend was also 
observed with UH-60A flight test correlation at // = 0.15 [9] 
and better correlation was obtained using a free wake 
geometry calculation method that combined the multiple- 
trailer wake with a simulation of the tip vortex formation 
process (consolidation). Calculations are carried out with the 
multiple-trailer wake with consolidation in order to identify 
whether the same improvement can be obtained at this low 
speed condition. Figure 22 shows section normal force and 
oscillatory flap bending moment correlation. The multiple 
trailer with consolidation model reduced the magnitude of the 
pulse due to BVI compared to the baseline rolled-up wake 
on the advancing side and, thus, improved the normal force 
correlation. A slight improvement is also observed on the 
front of the rotor disk. A similar improvement is obtained 
for the flap bending moment using the prediction with the 
multiple trailer and consolidation model. This prediction is 
very close to the coupled analysis. 

Structural loads calculation with prescribed measured 
airloads 

This section investigates the structural loads calculated from 
prescribed measured airloads (normal force, chord force, and 
pitching moment) using CAMRAD II alone. Good correlation 
between prediction and test depends on three factors: accurate 
airloads measurement, accurate structural loads measurement, 
and accurate structural model. Prescribing measured airlo ads 
permits validation of the structural model without error 
associated with aerodynamic modelings. This approach 
has been proved successful for the correlation with the 
NAS A/ Army UH-60A Airloads Program flight test data [19]. 
This study found the present structural model to be adequate. 

Structural loads are calculated for both /i = 0.3 and 0.4 and the 
results are compared with the measurements and the coupled 
analysis predictions as shown in Figs. 23 and 24. Good 
agreement between calculations and measurement is obtained 
for flap bending moment for both cases. Underprediction is 
observed for chord bending moment at /i = 0.3, but excellent 
correlation is obtained at /r = 0.4. Torsion moment correlation 
is poor. There are strong 5/rev responses, which were not 
observed with either CAMRAD II alone or coupled analysis. 
In order to understand poor correlation of torsion moments, 
section normal force and pitching moment near the blade 
tip area are compared in Figs. 25 and 26. Although there 
are pressure transducers installed at r/R = 0.965, airloads are 
not available at that radial location due to malfunction. One 
obvious outlier is r/R = 0.99 where there is a significant mean 
value change and a big hump in pitching moment in the first 
quadrant for both advance ratios. This may be the cause 
of a strong 5/rev response in the torsion moments. Without 
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measurements at r/R = 0.965, this is the only airload data 
available acting on the swept tip. Because of the sweep, these 
airloads will significantly twist the blade. The same issue was 
observed at other advance ratios (not shown in this paper) and 
thus assessment of structural loads with prescribed measured 
airloads was not successfully carried out in full. Detailed 
examination of pressures at this radial location is warranted. 

Parametric studies of maximum thrust condition 

At the high thrust condition, the predicted airloads featured 
large, first quadrant normal force and pitching moment 
variations that appear to stimulate a strong 5/rev structural 
response. These phenomena are not present in the measured 
data. In an effort to understand the cause of these 
mispredictions the effects of collective and pitch link stiffness 
on airloads and structural loads are investigated. All the 
results shown in Figs. 27 through 29 were obtained with the 
coupled analysis. Figure 27 shows the calculated pitch angle 
at the maximum thrust. Figure 27(a) is the result of two 
different collective angles, 11.88 deg and 14.85 deg. The 
11.88 deg was chosen to closely match with the measured 
collective and 14.85 deg was chosen because it is roughly 
the same collective change from the baseline collective of 
13.34 deg. The lower collective produced CtIo = 0.1187 
which is 5.42% lower than the measured value and the higher 
collective produced Ct/c = 0.1273 which is 1.43% higer. 
Larger collective required larger cyclic pitch angles to trim 
the rotor. Another variation, shown in Fig. 27(b), is pitch 
link stiffness with the baseline collective. The baseline pitch 
link stiffness value of 1090 ft-lb/deg, which is widely used for 
previous UH-60A analyses, was determined experimentally 
from direct measurements of the aircraft hardware by Kufeld 
and lohnson [20]. The soft value was obtained from the study 
by Shanley [21] under a NASA contract. 

Figure 28 shows section airloads due to the parametric 
changes. The lower collective (matching measured collective) 
reduces the normal force and pitching moment pulse in 
the first quadrant, but underpredicts negative normal force 
peak on the advancing side and misses the phase of two 
pitching moment peaks on the retreating side. The higher 
collective generates advancing blade stall spikes and degrades 
correlation. The soft pitch link improves correlation by 
reducing the normal force and pitching moment pulse in 
the first quadrant and correcting normal force phase on the 
advancing side. None of these improves pitching moment 
correlation across the front of the rotor disk. 

Figure 29 shows the effects of collective and pitch link 
stiffness on structural loads. Both lower collective and soft 
pitch link stiffness cases show slightly better correlation 
than the baseline results. The higher collective generates 
substantially larger 5/rev structural response. This is clearly 
illustrated by the 5/rev torsion moment harmonic correlation 
shown in Figs. 29(e) and 29(f). As the blade 2nd chord 


frequency is about 4.69/rev and first torsion frequency is about 
4.53/rev, the structural responses to the 5/rev airloads due to 
dynamic stall are magnified. In this case, the soft pitch link 
moves the first torsion frequency from 4.53/rev to 4.0/rev and 
prevents this excitation. 

Near the rotor thrust limit the measured thrust showed reduc ed 
sensitivity to collective change. The higher collective case 
overpredicted the rotor lift by just 1.43% but produced a 
dynamic stall event on the advancing side and significantly 
larger 5/rev structural responses. 

Conclusions 

Wind tunnel measurements of the rotor trim, blade airloads, 
and structural loads of a full-scale UH-60A Black Hawk 
main rotor are compared with calculations obtained using 
the comprehensive rotorcraft analysis CAMRAD II alone and 
a coupled CAMRAD II/OVERFLOW 2 analysis. A speed 
sweep at constant lift up to an advance ratio of 0.4 and a thrust 
sweep at constant speed into deep stall are investigated. 

From this study the following conclusions were obtained: 

Speed sweep 

1) The comprehensive analysis with lifting-line aerodynamics 
significantly overpredicted the magnitude of the down-up 
impulse in normal force on the advancing side at low advance 
ratio. However, the magnitude of the normal force is 
reasonably well predicted at higher advance ratios. Better 
normal force correlation is obtained at low advance ratio by 
using a free-wake geometry calculation method that combines 
the multiple-trailer wake with a simulation of the tip vortex 
formation process (consolidation). 

2) The coupled analysis shows significant improvement over 
comprehensive analysis. Pitching moment magnitudes are 
better captured in the coupled solutions. The shape of the 
airloads curves is better predicted with the coupled analysis. 

3) The comprehensive analysis shows reasonably good flap 
bending and torsion moment correlation but significantly 
underpredicts the magnitude of chord bending moment. The 
coupled analysis improves waveforms and phase of flap 
bending and torsion moment but does not improve the chord 
bending moment much. 

4) Structural loads were calculated by prescribing measured 
airloads. Strong 5/rev torsion moment, which were not 
observed in the test, were predicted and attributed to the lack 
of measurements at r/R = 0.965 and anomaly, especially in 
pitching moment, at r/R = 0.99. 

Thrust sweep 

5) The coupled analysis again shows significant improvement 
in the airloads predictions, especially capturing the magnitude 
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and phase of the two stall events in the fourth quadrant at the 
deeply stalled condition. 

6) Accurate analysis of the maximum thrust condition 
is challenging because large collective variation produces 
small thrust change. The higher collective case, although 
overpredicting lift by just 1.43%, produces a dynamic stall 
event on the advancing side and significantly larger 5/rev 
structural responses. 

7) As the blade 2nd chord frequency is about 4.69/rev and first 
torsion frequency is about 4.53/rev, the structural responses to 
the 5/rev airloads due to dynamic stall are magnified. In this 
case, a soft pitch link which moves the first torsion frequency 
from 4.53/rev to 4.0/rev improves the overall structural loads 
correlation. 
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Fig. 1 : UH-60A rotor system installed on the Large Rotor Test Apparatus in the NFAC 40- by 80-Foot Wind Tunnel. 
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(b) Strain gauge locations (□: flap bending, o: chord bending, and x: torsion 
moment) 


Fig. 2: Blade planform with locations of pressure transducers and strain gauges. 
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Fig. 3: Test conditions for speed and thrust sweep (circles represent test points investigated). 



(a) Near-body grid system (not every point shown) 



(b) Off-body grid system 


Fig. 4: OVERFLOW 2 grid system. 
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Fig. 7: Calculated and measured blade section normal force at Fig. 8: Calculated and measured blade section pitching 

r/R = 0.92 for speed sweep. moment at r/R = 0.92 for speed sweep. 
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Oscillatory FBM @ 40%R, ft-lb Oscillatory FBM @ 40%R, ft-lb Oscillatory FBM @ 40%R, ft-lb 








Fig. 9: Calculated and measured blade oscillatory flap 

bending moment at r/R = 0.4 for speed sweep. 


Fig. 10: Calculated and measured blade oscillatory chord 
bending moment at r/R = 0.4 for speed sweep. 
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Oscillatory TM @ 40%R, ft-lb Oscillatory TM @ 40%R, ft-lb Oscillatory TM @ 40%R, ft-lb 








Fig. 11: Calculated and measured blade oscillatory torsion 
moment at r/R = 0.4 for speed sweep. 


Fig. 12: Calculated and measured oscillatory pitch link load 
for speed sweep. 
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(a) 1/rev magnitude 



(d) 2/rev magnitude 



(g) 3/rev magnitude 
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Fig. 13: Comparison of harmonic magnitude of calculated and measured structural loads, fi = 0.3, C [J o = 0.09 
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Flap angle, deg Pitch angle, deg 



(a) y, = 0.3, C T /o- = 0.06 (b) /x = 0.3, C T /<r = 0.1255 

Fig. 14: Calculated and measured blade pitch angle for thrust sweep. 
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(a) fi = 0.3, C T /a = 0-06 


(b) fi = 0.3, C T /a = 0.1255 


Fig. 15: Calculated and measured blade flap angle for thrust sweep. 




Fig. 16: Calculated and measured blade section normal force at r/R = 0.92 for thrust sweep. 
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Fig. 17: Calculated and measured blade section pitching moment at r/R = 0.92 for thrust sweep. 




Fig. 18: Calculated and measured blade oscillatory flap bending moment at r/R = 0.4 for thrust sweep. 




(a) n = 0.3, C T /cr = 0.06 (b) n = 0.3, C T /a = 0.1255 


Fig. 19: Calculated and measured blade oscillatory chord bending moment at r/R = 0.4 for thrust sweep. 
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M C Oscillatory pitch link load, lb ►t) Oscillatory TM @ 40%R, ft-lb 




20: Calculated and measured blade oscillatory torsion moment at r/R = 0.4 for thrust sweep. 




Fig. 21: Calculated and measured oscillatory pitch link load for thrust sweep. 



(a) Blade section normal force at r/R = 0.92 



Fig. 22: Effects of wake modeling, /i = 0.2, Cl/ct = 0.92. 
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(b) Chord bending moment 


(b) Chord bending moment 




Fig. 23: Structural loads at 40%R calculated with prescribed 
measured airloads, n = 0.3, Cl/ct = 0.09 


Fig. 24: Structural loads at 40%R calculated with prescribed 
measured airloads, n = 0.4, Cl /<x = 0.09 
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Pitch angle, deg 




Fig. 25: Measured blade section airloads, /. l = 0.3, Cl/ct = 0.09 




Fig. 26: Measured blade section airloads, /x = 0.4, Cl/ct = 0.09 
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(a) Pitch angle with collective variation 


(b) Pitch angle with pitch link stiffness variation 


Fig. 27: Pitch angle and thrust, /x = 0.3, Ct / <J = 0.1255. 
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(a) Section normal force with collective variation (b) Section normal force with pitch link stiffness variation 




(d) Section pitching moment with pitch link stiffness 
variation 


Fig. 28: Blade section airloads atr/R = 0.92, /j = 0.3, Ct/o = 0.1255. 
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5/rev TM, ft-lb Oscillatory TM @ 40%R, ft-lb Oscillatory FBM @ 40%R, ft-lb 








Fig. 29: Oscillatory structural loads, /r = 0.3, Ct/ct = 0.1255. 
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